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Abstract 

The Bak-Chen-Tang forest fire model (1) was proposed as a toy model of turbulent 
systems, where energy (in the form of trees) is injected uniformly and globally, but 
is dissipated (burns) locally. We review our previous results on the model (2; 3) and 
present our new results on the statistics of the higher-order moments for the spatial 
distribution of fires. We show numerically that the spatial distribution of dissipation 
can be described by Log-Poisson statistics which leads to extended self-similarity 
(ESS) (4; 5). Similar behavior is also found in models based on directed percolation; 
this suggests that the concept of Log-Poisson statistics of (appropriately normalized) 
variables can be used to describe scaling not only in turbulence but also in a wide 
range of driven dissipative systems. 
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1 Introduction 



Understanding the spatial distribution of energy dissipation is crucial in the 
study of homogenous, isotropic, fluid turbulence. Kolmogorov, in his cele- 
brated 1941 theory (K41) (p), first conjectured that, in the inertial range, 
the velocity structure functions scale as a power law in the separation I. The 
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theory implicitly assumes uniform energy dissipation characterized by a mean 
dissipation rate. Experiments from the recent decades showed clear devia- 
tions from the K41 scaling, which are believed to be related to the highly 
intermittent and non- uniform spatial distribution of dissipation. JIB Man^ 
phenomenological models of energy cascade and intermittency have been pro- 
posed to understand the deviations from the K41 sc aling . The most famous are 
the Lognormal model (0), the Multifractal model (flpf). and the Log-Poisson 
model by She and Leveque (SL) (fill). The SL model, based on a hypothesis 
about the behavior of the moments of the energy dissipation, leads to a predic- 
tion of the exponents for the velocity structure functions that appears to be in 
excellent agreement with experimental results. Benzi et al. (0) discovered an 
important property they termed Extended Self-Similarity (ESS): the p th order 
structure function has a power-law dependence on the third-order structure 
function. They found that the validity of ESS extended almost down to the 
dissipative range, roughly five times the Kolmogorov scale. The connection 
between the dissipation distribution and the scaling in the inertial range is 
further supported by these investigations. 

It remains a great challenge to understand the dynamical mechanism for the 
inertial range scaling and ESS in turbulence. So far it has not been feasible 
to attack the problem directly from the underlying equations for turbulence. 
Instead we choose to study a few simple models of "turbulent systems" , such 
as the BCT forest fire model. We hope that, by illustrating the emergence of 
complex spatial distribution of dissipation and ESS in such simple dynamical 
models, we can gain some intuition about the dynamical mechanism for energy 
cascade and scaling in fully developed turbulence. 



2 Scale-Dependent Dimension in the BCT Forest Fire Model 

The BCT forest fire model is defined as follows. On a d- dimensional lattice of 
size L d , trees, representing energy, are grown randomly at empty sites at a rate 
p. At each time step, trees that are on fire are burnt down (the site becomes 
empty at the next time step) and neighboring trees are ignited. The fires die 
out when there are no more trees in their neighborhood. These processes repre- 
sent the dynamics of energy input and dissipation. Numerical studies indicate 
that there is a correlation length £ (£ ~ p~ 2 / 3 in 3D) in the system (J2J). If L < £ 
the fires cannot be sustained, so we consider only the case L > £. Larger and 
larger systems can be studied in the limit p — > 0. Despite the simplicity of the 
model, many fascinating scaling behaviors emerge. For example, Broker and 
Grassberger (flih found anomalous scaling in three and four dimensions; in 
particular, they found periodic oscillations in autocorrelation functions with 
a period T which diverges as 1/V 3 with (3 ~ 0.77 in three dimensions. 
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Fig. 1. D(l) vs log 2 I for the BCT forest fire model with L = 1024 and p = 0.0003. 
Averages were obtained using 10 6 time steps (after skipping 200,000 time steps). 
The fits using the expression of scale-dependent dimension and the simple "fractal" 
scaling with crossover to 3D are also shown. 

One of the interesting features of the model is the unconventional scaling 
that can be interpreted geometrically as a scale-dependent fractal dimension: 
the apparent dimension D(l) characterizing the distribution of fires at the 
length scale I in three spatial dimensions increases linearly in log(7) (0; 0). 
The dimension increases gradually as one steps further and further backwards, 
and views the system at larger and larger length scales. At some point, the 
correlation length is reached; the distribution of fires remains uniform beyond 
that length. This is the main geometrical picture for the spatial distribution of 
the fires. We have repeated the simulation for larger lattices (up to 2048 3 ) and 
smaller p (as small as 0.00015); so far no deviation from the scaling described 
above is found. As was done in the previous simulation (0), we determine the 
mean number of fires n(l) within boxes of size I < £ that contain fires. The 
linear dependence of D{1) on log(Z) corresponds to 





which gives rise to 




(2) 



3 



20 




2 4 6 8 10 



log 2 l 

Fig. 2. log 2 eoo(/) vs log 2 1 for the BCT forest fire model with L = 1024 and 
p = 0.0003 

Fig. 1 shows the numerical values of D(l) vs log 2 /. Here D(l) is calculated using 
the numerical derivative: D(2 1 ^ 2 l) = log 2 (n(2Z)) — log 2 (n(7)). In the figure we 
also show the two-parameter fit using the above equation. There is excellent 
agreement between the data and the our proposed expression of D(l) (even a 
one-parameter fitting with Iq fixed to Iq — 1 is very good). For comparison we 
also show the fitting using a simple "fractal" scaling: n(l) = al b + cl 3 which 
gives D(l) = (abl b + 3c/ 3 )/(a/ b + cl 3 ). This three-parameter fitting is very- 
bad. It is interesting to note that the distribution of luminous matter in the 
universe might also follow a scaling form similar to the one given in Eq. 2 (0). 



3 Extended Self-Similarity in the Forest Fire Model 



The scale dependent dimension shown in the previous section only concerns the 
first moment of the fire distribution. We will show that the higher-order mo- 



ments also exhibit interesting scaling behaviors. Preliminary study in Ref. (|14h 
indicated that the higher-order moments can be described by extended self- 
similarity used to describe the moments of energy dissipation in turbulence. 
We present here a detailed study of the higher-order moments. 

We start with the higher-order moments given by 

M q {l) =< m{l) q >, (3) 
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where m(l) is the number of fires contained in a box of linear size I] we empha- 
size that the average is again over all the boxes that contain fires. Following 
Dubrulle's analysis of generalized scale covariance and ESS in turbulence (|5|), 
we consider the normalized variables, defined by tt(£) = m(/)/e 00 (/), where 
eoo(Z) given by 

M q+1 (l) 

eoo(Z) = hm . (4) 

g^oo M 9 (/) 



We focus on these normalized variables and calculate their moments, defined 
as 

n q (i) =< > . (5) 



We will show numerically that the scaling behavior of {IT^} is consistent with 
the Log-Poisson statistics. Let it (I) = fl 3 . If j is a random variable described by 

a Poisson distribution Pi{j) = — , then 7r(7) is described by a Log-Poisson 
distribution. The moments of 7r(Z) can then be written as 

n 9 (0 = E^'^r^ = exp w - 1)]. (6) 

j=Q J' 



It is easy to see that the moments from this Log-Poisson distribution are 
related by 

n,(0 = n p (/)«^ = n 1 (O f(<?) , (7) 



where 

C0?) = (l-/3 9 )/(l-/3). (8) 

Note that ((q) is not linear in q as would be the case for simple scaling. The 
moments of the original variables m(l) are given by 

M a (0 = 4 > (0ni(0 f(9) . (9) 



To test whether the moments can be described by the Log-Poisson distribution 
numerically, we first fit the data M q (l) to the above expression. The actual 
fitting is done using the following expression: 

log 2 (M ?+1 (Z)/M ? (0) = o(Z) + b(l)(3\ (10) 
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Fig. 3. Uq vs II3. The lines are drawn with slopes Cq/Cz-, where Cq are defined in 
Eq. 8. The value of (3 = 0.75 is chosen as explained in the text. 

where a(l) = log 2 e^Q) and b(l) = log 2 (Ili(/)). We first determine (3 as follows: 
For each value of I we do a linear fitting of log 2 (M (?+ i(/)/M IJ (7)) vs (3 q , and 
calculate the overall error, which depends on (3. The best value of /3, corre- 
sponding to the minimum overall error, is found to be (3 ~ 0.75. With [3 fixed 
to this value, we find a(l) or e oc (Z) from the fitting; this is plotted in Fig. 2. 

For I < £ (£ is about 2 8 for the value of the parameter p = 0.0003 used in the 
simulation), e^l) ~ Z 1-8 . Note that if the original variable m(l) is defined as 
the density of fires ( instead of the number of fires) in a box of linear size /, e 00 (Z) 
will be approximately given by l 18 l~ 3 = Z~ L2 — but the normalized variables 
remain the same. Given e^/) we can evaluate the normalized moments H q (l)- 
If the distribution is Log-Poisson, then Tl q vs Tl p in the log-log plot should fall 
on a straight line with a slope given by C(?)/C(p)- Fig- 3 shows U g vs il 3 for 
9 = 1,2,- -.,10. 

Indeed, we find straight lines in the log-log plots with slopes matching ((q)/((3). 
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Fig. 4. D(l) vs log 2 I for 2 + 1 DP model with L = 8192 and p = 0.28734. 5 x 10 5 
time steps (after skipping 100, 000 time steps) were used to obtain the averages. 

This is extended self-similarity (in the normalized variables) for the forest fire 
model. This is particularly interesting, given the fact that there is actually no 
scaling regime in the plot of M g (l) or U q (l) vs I. The Log-Poisson distribution 
provides an excellent description of the normalized moments. 



4 Extended Self-Similarity in Models of Directed Percolation 

It is of interest to determine whether the two properties of the forest fire 
model described in the previous section are displayed by other models. The 
first, the special form of the scale-dependent dimension appears to be unique 
to the forest fire model; so far we have not found other types of models which 
exhibit similar scaling in the first moment. However, the second property of 
extended self-similarity in normalized variables as described in the last section 
appears to be much more general. Since the fire propagates and spreads in 
the background of trees which has presumably a hierarchical and certainly 
a dynamically changing topology it is interesting to consider simpler models 
which describe a spreading process. The most obvious one is the well-studied 
model of directed percolation (DP), (flih We found to our surprise that quite 
a few models of directed percolation (DP) and their variants exhibit a similar 
form of ESS based on log-Poisson statistics. As an example we consider 2 + 1 
DP in the supercritical phase p > p c , which is very different from the 3d forest 
fire model. Fig. 4 shows D(l) vs I for this model. Here we use a lattice of size 
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Fig. 5. U q vs II3 for 2 + 1 DP model. The lines are drawn with slopes Cg/C3> where 
C q are defined in Eq. 8. The value of (3 = 0.76 is used. 



8192 x 8192 and p = 0.28734. Over half millions time steps were used to obtain 
the averages. 

There is an approximately scale-invariant regime corresponding to D(l) being 
constant, in contrast to the 3d forest fire model. However, the normalized 
higher order moments can still be described the Log-Poisson statistics (with 
( q defined in Eq. 8) and the associated ESS, as can be seen from Fig. 5. 
This data provides evidence that either the model studied or a model in its 
vicinity in parameter space exhibits the behavior displayed asymptotically; the 
proximity to the appropriate fixed point controls this behavior and elucidating 
its features would be of interest. 



5 Summary 



We presented extensive numerical evidence to show that the BCT forest fire 
model exhibits two interesting properties: a logarithmic scale-dependent di- 
mension and Log-Poisson distribution of normalized variables (with ( q not 
linear in q) and associated extended self-similarity. The former is likely a 
unique feature of the forest fire and closely related models; the latter, however, 
appears to be more general. We have shown that models based on directed 
percolation also exhibit Log-Poisson statistics (for appropriately normalized 
variables) similar to that in the forest fire model. This suggests that the Log- 
Poisson distribution and the associated ESS can be used to describe scaling 
in a wide range of driven dissipative systems. We hope that studying sim- 
ple models such as the forest fire model will provide us with some intuition 
and physical pictures that might eventually be useful for understanding the 
statistics of energy dissipation and scaling in fully developed turbulence. 

This work was supported by the National University of Singapore research 
grant R-151-000-028-112. KC thanks Maya Paczuski, Peter Grassberger, and 
Itamar Procaccia for helpful discussions. This work was originally initiated by 
Per Bak, who was very interested in scaling and ESS in turbulence. He had 
hoped that a physical picture for ESS could emerge in a simple dynamical 
model such as the forest fire model. It is sad that he is not around to offer 
his unique insight on these topics. For many, particularly those who worked 
closely with him, his premature death has left a big void. 
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